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ABSTRACT 

A comprehensive census of Galactic open cluster properties places unique constraints 
on the Galactic disc structure and evolution. In this framework we investigate the 
evolutionary status of three poorly-studied open clusters, Berkeley 31, Berkeley 23 
and King 8, all located in the Galactic anti-centre direction. To this aim, we make use 
of deep LBT observations, reaching more than 6 mag below the main sequence Turn- 
Off. To determine the cluster parameters, namely age, metallicity, distance, reddening 
and binary fraction, we compare the observational colour-magnitude diagrams (CMDs) 
with a library of synthetic CMDs generated with different evolutionary sets (Padova, 
FRANEC and FST) and metallicities. We find that Berkeley 31 is relatively old, with 
an age between 2.3 and 2.9 Gyr, and rather high above the Galactic plane, at about 
700 pc. Berkeley 23 and King 8 are younger, with best fitting ages in the range 1.1-1.3 
Gyr and 0.8-1.3 Gyr, respectively. The position above the Galactic plane is about 500- 
600 pc for the former, and 200 pc for the latter. Although a spectroscopic confirmation 
is needed, our analysis suggests a sub-solar metallicity for all three clusters. 

Key words: Hertzsprung-Russell and colour-magnitude diagrams, Galaxy: disc, open 
clusters and associations: general, open clusters and associations: individual: Berkeley 
31, open clusters and associations: individual: Berkeley 23, open clusters and associa- 
tions: individual: King 8 



1 INTRODUCTION 

This paper is part of the BOCCE (Bologna Open Clus- 
ter Chemical Ev o lution ) project, described in detail by 
iBragaglia fc Tosil |2006) and aimed at deriving precise 
and homogeneous ages, distances, reddenings, and chemi- 
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tuto Nazionalc di Astrofisica, Italy; LBT Beteiligungsgcscllschaft, 
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cal abundances for a large sample of open clusters (OCs). 
Our final goal is to study the disc of our Galaxy and 
its formation and evolution. In fact OCs are among the 



best t racers of th e disc properties (e.g., Panagia fc Tosil 
I l98l| : lErieil Il995l. iTwarog. Ashman, k, Anthonv-Twaroel 
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1 19971 : iFreeman fc Bland-Hawthornll2002l ). We have already 
published results based on photometry for 23 OCs (see 
IBragaglia fc Tosil2006l ; lAndreuzzi etld1l201ll . and references 
therein), concentrating on the old ones, the most important 
to study the early epochs of the Galactic disc. With the 
three OCs presented here we have a sample of 20 clusters 
with ages older than 1 Gyr, i.e. about 10 per cent of all 
known old clusters (see the catalogue bv lDias et al. 2002b 
and its web updates). 

The three clusters examined in this paper are Berke- 
ley 23 (I = 192.6°, b = 5.4°), Berkeley 31 (I = 206.2°, b = 
5.1°), and King 8 (I = 176.4°, b = 3.1°). They are old, dis- 
tant clusters in the anti-centre direction. They were selected 
because they should all lie beyond a Galactocentric distance 
of 12 kpc and could then be useful to understand the nature 
and properties of the outer Galactic disc. In particular, they 
are located in the region where the radial metallicity distri- 
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bution seems to change its slope (e.g.jFriel et"afll201Cl : sec 
the discussion in I Andreuzzi et al]|201ll ). 

These three clusters have already been studied to dif- 
ferent degrees in the past, but with contrasting results. We 
present here high quality photometric data, to improve on 
the previous determinations of their parameters. As done for 
all our past work we use comparison of observational CMDs 
to synthetic ones generated using different sets of evolution- 
ary tracks. 

B erkeley 23 h as be en the subject of two photometric 
works. lAnn et a l. (2002) observed it with a 1.8m telescope 
in UBVI as part of a study of 12 OCs, deriving a reddening 
E(B - V) = 0.40 ± 0.05 (from the two-colour diagram), 
a distance modulus (m — M)o = 14.2 ± 0.3, a metallic- 
ity [Fe/H]=+0.07, and an age of 0.79 Gyr. However, the 
isochro ne they chose does no t seem to reproduce the red 
clump. Hasegawa et al. I (|2004 ) observed 14 OCs with a 65cm 
telescope, obtained BVI photometry and derived the clus- 
ters' parameters using isochrones. Their results for Be 23 
differ from the previous ones: E(B — V) = 0.30, (m — M)o = 
13.81, metallicity Z=0.004 (equivalent t o [Fe/ Hl~ -0.7), 
and age 1.8 Gyr. lAhumada fe Lapassetl l|2007h . who pre- 
sented a catalogue of Blue Straggler stars in OCs (they have 
about 1900 candidate BSS in about 430 clusters), find that 
Be 23 is rich in this kind of stars. 

Ber k eley 31 is the best studied of the three OCs. 
iGuetteri (|l993T ) published UBVI photometry obtained at 
the lm USNO telescope. He determined the reddening 
and the metallicity using the two-colour diagram, and 
the distance and age with isochrone fit, deriving E(B — 
V) = 0.13, (m - M) = 13.6, [Fe/H]= -0.4, and age 
8 Gyr. This old age is probably an artifact of the out- 
dated stellar models. Furthermore, the isochrones shown 
do not reach the red clump, while the difference in mag- 
nitude between it and the main sequence (MS) turn- 
off (TO) is a powerful age indicator (for different defini- 
tions of this SV see, e.g., ITwarog fc Anthony- Twarodl 1989 ' 



ell 19891 ; 
(|l994T ) 



iPhelps. Janes, fc Montgomery! I1994T ). iPhelps et all (|l994 
obtained VI photometry, but could not derive directly 8V; 
the value they give would put Be 31 at the same age of 
M67. [Hasegaw a~et al.l (|2004T ) ■ who call it Biurakan 7, derive 
from isochrone fit E(B - V) = 0.15, (m - M) = 14.83, 
Z=0.008, and age of 2.2 Gyr. The discording results are 
not limited to the photometry. iFriel et all l|2002ft . using low 
resolution spectra of 24 stars (17 defined as members) de- 
rived an average radial velocity < RV >= +41 (rms=15) 
km s -1 and a metallicity [Fe/H] =— 0.40 ( rms 0.16) dex. 
lYong. Carney, fc Teixera de Almeidal (|2005h obtained high 
resolution spectra for five stars; they notice that the RV dis- 
persion is very high for an OC, maybe because they observed 
by chance a large fraction of binaries or because their sam- 
ple includes field stars with RV similar to the one of Be 31. 
They were able to perform an abundance analysis only on 
one star a nd derived [ Fe/H]= —0.57. A similar situation was 
found by IFriel et all (|2010l l who obtained high resolution 
spectra of another two st ars, with discrepan t RVs -differing 
also from the values by lYong et al.l ([2005)- and metallic- 
ity (although the two values, [Fe/H] = -0.22 and = -0.32 
are not too distant). They could not determine which of the 
two stars, if any, is a true cluster me mber. They also noticed 
that t he published proper motions (|Dias. Lepine. fc Alessil 
2002a) should be of field stars in the direction of the clus- 



ter, given their bright magn itudes. Finally, also for Be 31 
lAhumada fc Lapassetl (|2007r i found a large number of BSS. 

The first photometry ( i n the GRU filters) of King 8 was 
published by ISvolopoulosI (|l965T ). who described a young 
cluster, with an age of a bout 5 Myr, E(B — V) = 0.34, 
and (m - M) = 13.3. IChristianl (|l98ll ) presented BV 
photographic photometry and low resolution spectra. She 
gave E(B — V) — 0.7 (with a possible differential red- 
dening up to 0.1 mag), Z< 0.001, distance from the Sun 
of 3-4 kpc, and age 0.8 Gy r. These values were partially 
changed by IChristianl (|l984r ), where E(B - V) = 0.55 and 
[Fe/H]~ —0.4 or —0.5 dex. These two works were later 
used by ITwarog et al.l |T997) to derive a larger distance, 
( m — M) = 15.3 0, adop ting the metallicity [Fe/H]= -0.46 
of IFriel fc Janes (1 19931). based on two stars observed at 
low resolution. Koposov. Glushkova. fc Zolotukhml (120081 ) 
re-discovered King 8 in their automated search for new 
OCs on the 2MASS images. They derived the cluster pa- 
rameters using isochrones in the J, J — H plane, finding 
E(B — V) = 0.44, distance from the Sun of 2.4 kpc, and 
age 1.1 Gyr. They al so noted that their numbers differ from 
what is listed in the lDias et a l. (2002b) catalogue. 

In summary, even if apparently these three OCs have 
already been studied, their properties are stil l insecure. 
Be 2 3 could be rather young and metal rich (|Ann et al.l 
120021 : 0.79 Gyr, [Fe/H1=+0 .07) or rather old and metal- 
poor (jHasegawa et al.ll2004l : 1.8 Gyr, [Fe/H]=-0.7). Be 31 
could be very old or of inte rmediate age, and at very dif- 
ferent Galactocentric radii l|Guetten[l993l : 8 Gyr, 13 kpc; 
lHasegawa et al.ll2004l : 2.2 Gyr, 17 kpc), while there is some 
concordance on E(B — V) around 0.13-0.15 and metallicity 
around —0.4 dex (but with caveats on membership). King 8 
has probably an age around 1 Gyr and sub-so lar metallicity 
but t he reddening and distance are un certain ( Twaro g et al.l 
1 19971 : E(B - V)=0.55, d Q =5.24 kpc: iKoposov et al.ll200&!T 
E(B - V)=0.U, d =2.4 kpc). 

In the present paper we will describe our new observa- 
tions and the resulting CMDs (Sect. 2), the radial extension 
of the clusters (Sec. 3), the derivation of their age, distance, 
reddening, and metallicity using comparison to synthetic 
CMDs (Sect. 4). A discussion and summary is presented 
in Sec. 5. 



2 THE DATA 
2.1 Observations 

The three clusters were observed in service mode at the 
Large Binocular Telescope (LBT) on Mt. Graham (Arizona) 
with the Large Binocular Camera (LBC) in 2008. There are 
two LBCs, one optimised for the UV-blue filters and one 
for the red-IR ones, mounted at each prime focus of the 
LBT. Each LBC uses four EEV chips (2048x4608 pixels) 
placed three in a row, and the fourth rotated 90deg with 
respect to the others. The field of view of LBC is equiva- 
lent to 23'x23', with a pixel sampling of 0.23". For tech- 
nical reasons, only LBC-Blue was available at the time of 
our observations, so we collected only B and V data. The 
clusters were positioned in the central chip (# 2) of the 
LBC-Blue CCD mosaic. Table [T] gives a log of the observa- 
tions. While the seeing was good (below 1" for all images, 
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Table 1. Log of observations for the LBT and TNG, respectively. The exposure times are in seconds. 



Cluster 


RA 


Dec 


Date 


B 


V 




(2000) 


(2000) 




exptime 


exptime 








LBC@LBT 






Be 23 


06 33 15 


+20 34 50 


02 Dec 2008 


Is, 3x5s, 4x90s 


Is, 3x5s, 3x60s 


Be 31 


06 57 37 


+08 21 18 


02 Dec 2008 


Is, 3x5s, 3x90s 


Is, 3x5s, 3x60s 


King 8 


05 49 24 


+33 40 60 


02 Dec 2008 


Is, 3x5s, 4x90s 


Is, 3x5s, 3x60s 








DOLORESOTNG 






Be 23 


06 33 15 


+20 31 30 


03 Jan 2009 


10s, 20s 


5s, 10s, 20s 


Be 31 


06 57 37 


+08 17 20 


03 Jan 2009 


10s, 20s 


10s, 20s 


King 8 


05 49 18 


+33 37 50 


03 Jan 2009 


10s, 20s 


10s, 20s 



and below 0.8" for many), the observing conditions were not 
photometric. To calibrate our photometry we obtained a few 
shallow images using DOLORES (Device Optimised for for 
the LOw RESolution) at the Italian Telescopio Nazionale 
Galileo (TNG) on a night dedicated to another programme 
and to service observations. For this instrument the field 
of view is 8.6' x 8.6', with a scale of 0.252 "/pix. The see- 
ing was bad (about 3") but this is not a problem since the 
fields are not crowded. Only the central part of the clusters 
was observed, barely reaching stars at the turn-off of the 
main sequence. The t wo standard st ar fields PG0918+029 
and PG1323-085 from lLandoltl l| 19921 ) were also obtained. 

2.2 Data reduction 

Given the large area covered by LBC and the small angular 
dimensions of our targets (see Sections [3]|4| we only reduced 
two of the four chips, the one centered on the cluster (chip 
#2) and a second one to be used as a comparison to separate 
the cluster from the field stars (chip #1). 

The raw LBC images were corrected for bias and flat 
field, and the overscan region was trimmed using a pipeline 
specifically developed for LBC image prereduction by the 
Large Survey Center team at the Rome Astronomical Obser- 
vatorjQ. The source detection and relative photometry was 
performed independently on each B and V i mage, using the 
PSF- fitting code DAOPHOTII/ALLSTAR ljStetsonlll987l . 
1 19941 ). The brightest stars, saturated in the deepest images, 
where efficiently recovered through the accurate magnitude 
measures from the short exposure time images. The average 
and the standard error of the mean of the independent mea- 
sures obtained from the different images were adopted as the 
final values of the instrumental magnitude and uncertainty. 

More than 300 stars from the Guide Star Catalogue 2.3 
where used as astrometric standards to find an accurate as- 
trometric solution to transform the instrumental positions, 
in pixels, into J2000 celestial coordinates. To this aim wc 
adopted the code CataXcorr, developed by Paolo Monte- 
griffo at INAF - Osservatorio Astronomico di Bologna, and 
successfully used by our group in the past 10 years. The 
r.m.s. scatter of the solution was ~ 0.3" in both RA and 
Dec. 

We derived the completeness level of the photometry 
by means of extensive artificial stars ex p erime nts following 
the recipe described in iBellazzini et alJ l|2002t ). More than 



150000 artificial stars have been uniformly distributed on 
the chip sampling the cluster in groups of 9000 stars at a 
time, to avoid changing the crowding conditions. The whole 
data reduction process has been repeated as in the real case 
and the fraction of recovered stars was estimated at each 
magnitude level. We also checked that the completeness does 
not change from the central to the external parts. Results 
for the three clusters are presented in Table [2] 

With the standard stars observed at the TNG we de- 
rived the following calibration equations 

B = b + 0.1024 x (6 - v) + 1.3580 (rms = 0.0118) 

V = v - 0.0506 x (6 - v) + 1.0416 (rms = 0.0167) 

where B, V are the calibrated magnitudes and b, v the in- 
strumental ones (after aperture correction). Finally, the deep 
LBT catalogues were cross-correlated with the calibrated, 
shallow TNG ones and the stars in common were used 
as "secondary standards" to extend the calibration to the 
whole catalogues, both for the central and the external chip. 

We produced three catalogues with identification, equa- 
torial coordinates, B and V magnitudes, with errors contain- 
ing about 3600, 6500, and 4000 stars for Be 23, Be 31, and 
King 8, respectively. These catalogues will be made public 
through the WEBDA0. 

Fig. [T] shows the resulting CMDs from our photometry 
for the central and comparison field chips (upper and lower 
panels, respectively). In the upper panels we also indicate 
the photometric errors in magnitude and colour as derived 
from extensive artificial star tests. They range from less than 
0.01 mag at the bright limit to less than 0.1 mag around 
V = 24. The three OCs are visible over the field background, 
but they represent a minority of the stars in each frame. In 
the most favourable case, Be 31, the main sequence is rather 
tight and clearly extended down to V « 23 — 24, while for 
Be 23 and King 8 the MS is recognisable from the field with 
some difficulty. A discussion of the field population sampled 
in these CMDs is presented in Sect. 5. 

2.3 Comparison to previous data 

We comp a re ou r photome try for t he t hree OCs with 
Ann et al.l (I2002T ) for Be 23. iGuetterl (|l993l ) for Be 31, and 
ChristianTl|l98ll ) for King 8. We do not c onsider the pho- 
tometry of Be 31 by iPhelps et al. I (1 19941 ) because it has 



1 http://lsc.oa-roma.inaf.it/ 



2 http:/ /www. univie.ac.at/webda/webda. html 
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Figure 1. Upper panels: CMDs of Be 23, Be 31, and King 
Photometric errors, as derived from the artificial star tests. 



3 (chip #2). Lower panels: CMDs of the relative comparison fields (chip #1). 
are shown on the right side of the CMDs. 



only V and / and that of lHasegawa et all (|2004l ) of Be 31 
and Be 23 because they made public through the WEBDA 
only the bright stars (V 18). 

We downloaded the literature photometry files from the 
WEBDA site and cross-correlated our catalogues with them. 
Figs [U121 and 2] show the results of this comparison between 
our and literature magnitudes. 

For Be 23 (Fig. [2| the difference in the B filter are 
very small (on average -0.013 mag), while the V shows a 
small trend with magnitude. We do not know if this is due 
to our photometry or theirs; however, no trend is present 
for the two other OCs, that were calibrated with the same 
equations, and this is in favour of the latter possib ility. In the 
same figure we also show the CMD obtained bv lAnn et al.l 



(2002) to give an immediate impression of the much better 
quality of our data (see Fig. [1} . 

In the case of Be 31 (Fig. [3} , we find small average differ- 
ences both in B and V (0.02 9 and 0.058 ma g, respectively), 
and no trends. The CMD bv lGuetterl ljl993h . even if of good 
quality, is poorer than ours, which reaches fainter and shows 
details (e.g. near the MS TO) not visible in Guetter's. 

Finally, King 8 (Fig. [4j has only the photographic data 
bv lChristianl (|l98ll ). While the difference between the two 
B values are small (-0.030 mag), for V we reach an average 
-0.185 mag. Again, this cluster was calibrated exactly as 
the other two and, given the much smaller differences with 
completely independent sources found for the other OCs, we 
think that the problem is (mostly) in the old photographic 
photometry. 
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Table 2. Completeness of the photometry for the three clusters. 



mag 


compl B 


compl V 


compl B 


compl V 


compl B 


compl V 








Be 23 






Be 31 








King 


8 




17.25 


1.00 


±0.11 


0.99 


±0.05 


0.95 


±0.11 


0.96 


±0.06 


0.98 


±0.05 


0.99 


±0.02 


17.75 


1.00 


±0.11 


0.98 


±0.04 


0.95 


±0.12 


0.97 


±0.06 


0.98 


±0.05 


0.99 


±0.04 


18.25 


0.99 


±0.06 


0.98 


±0.04 


0.98 


±0.09 


0.98 


±0.05 


0.99 


±0.04 


0.99 


±0.03 


lo. / 


0.98 


±0.04 


0.97 


±0.03 


0.97 


±0.06 


0.97 


±0.05 


0.99 


±0.04 


0.98 


±0.03 


19.25 


0.97 


±0.04 


0.97 


±0.03 


0.98 


±0.05 


0.95 


±0.04 


0.98 


±0.03 


0.97 


±0.03 


19.75 


0.98 


±0.03 


0.97 


±0.03 


0.97 


±0.05 


0.96 


±0.04 


0.98 


±0.03 


0.97 


±0.03 


20.25 


0.97 


±0.03 


0.96 


±0.03 


0.96 


±0.04 


0.95 


±0.04 


0.98 


±0.03 


0.97 


±0.03 


20.75 


0.97 


±0.03 


0.95 


±0.03 


0.96 


±0.04 


0.94 


±0.04 


0.98 


±0.03 


0.96 


±0.02 


21.25 


0.98 


±0.03 


0.94 


±0.03 


0.95 


±0.04 


0.92 


±0.04 


0.97 


±0.03 


0.95 


±0.02 


21.75 


0.96 


±0.03 


0.93 


±0.03 


0.94 


±0.04 


0.90 


±0.03 


0.96 


±0.03 


0.95 


±0.02 


22.25 


0.95 


±0.03 


0.90 


±0.02 


0.94 


±0.04 


0.89 


±0.03 


0.96 


±0.02 


0.92 


±0.02 


22.75 


0.94 


±0.03 


0.87 


±0.02 


0.90 


±0.03 


0.84 


±0.03 


0.95 


±0.02 


0.90 


±0.02 


23.25 


0.91 


±0.02 


0.78 


±0.02 


0.88 


±0.03 


0.80 


±0.03 


0.93 


±0.02 


0.88 


±0.02 


23.75 


0.89 


±0.02 


0.55 


±0.02 


0.85 


±0.03 


0.72 


±0.02 


0.93 


±0.02 


0.80 


±0.02 


24.25 


0.85 


±0.02 


0.29 


±0.01 


0.82 


±0.03 


0.53 


±0.02 


0.89 


±0.02 


0.59 


±0.01 


24.75 


0.74 


±0.02 


0.10 


±0.01 


0.73 


±0.02 


0.26 


±0.01 


0.83 


±0.02 


0.26 


±0.01 
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Figure 2. CMD of Be 23 by Ann et al. (2002) on the same 
scale used for next figures with our data. Right panels: differences 
between our magnitudes and theirs, in B (upper panel) and V 
(lower panel). Open symbols indicate all stars in common, light 
blue filled ones the stars within 2a from the average, used to 
compute the mean differences. 
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Figure 3. Left panel: CMD of Be 31 by Guetter (1993) and 
comparison with our data (see previous figure) 
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Figure 4. Left panel: CMD of King 8 by Christian (1981) and 
comparison with our data (see previous figures). 



3 THE RADIAL PROFILES 

Taking advantage of the combination of the wide field and 
deep imaging capabilities of LBC on LBT, we determined 
the cluster structural parameters from the star density pro- 
files. As a first step, the centre of gravity (C gra v) was esti- 
mated simply by averaging the a and 8 coordinates of clus- 
ters stars. In order to avoid completeness effects and strong 
field contamination, we considered samples with different 
limiting magnitudes (V < 20, 20.5, 21) and with colour 
B — V < 1.5. For each sampl e we iteratively cal culated 
the barycentre of the stars (see iBeccari et al.l l200fjt ) inside 
a circular area of 70"radius. The size of the area, somehow 
arbitrary, allows us to sample a gradient in the stellar radial 
distribution (see details later on this paragraph) and to min- 
imise field contamination when not removed after colour and 
magnitude cuts. The centre of gravity has been calculated as 
the average of the measures in the three magnitude ranges 
and the standard deviation as the related error. The values 
of Cgrav and relative errors are reported on Table where 
we also indicate the coordinates found in the WEBDA, since 
they are slightly different. 

The projected density profile was determined for each 
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Figure 5. Observed surface density profiles (filled circles and 
error bars) in units of number of stars per square arc-second. The 
solid lines are the King models that best fit the observed density 
profile over the entire cluster extension. 



cluster usi ng direct star counts following the procedure de- 
scribed in lLanzoni et all (|201ttl . Shortly, the photometric 
catalogue is divided in concentric annuli centered on the 
cluster Cgrav Each annulus is then divided in at least three 
angular sectors. The density of the annulus is calculated as 
the mean of the stellar density in each sector with the stan- 
dard deviation of the mean used as the associated error. 
The observed radial density profile for each target cluster 
is shown in Fig. [5] with the abscissas corresponding to the 
mid-point of each radial bin. Notice that, in order to min- 
imise field contamination, we used only stars with magni- 
tudes V < 21 and colour B - V < 1.5. 

The cluster structural parameters can be estimated 
through a best fittin g procedure of the derived density pro- 
file with a suitable iKind (|l966l ) model (as done, e.g., in 
ISollima et al.ll2010l ). An isotropic, single-mass King model 
pro jected onto the cluste r area was computed adopting 
the ISigurdsson fc Phinnevl l|l995l ) code. The model produc- 
ing the best fit for each cluster is shown in Fig. [5] together 
with the structural parameters core radius (r c ), half-mass ra- 
dius {th), tidal radius (n) and concentration (c = \n(r t /r c ), 
all indicated in Table [3] The residual of the fit of the model 
with each observed point is shown in the lower panel of each 
density plot. Notice that the estimate of the background 
level is based on the star counts in the peripheral regions 
where the number of field stars is significantly higher than 
that of cluster stars. 



4 CLUSTERS' PARAMETERS VIA 

SYNTHETIC COLOUR-MAGNITUDE 
DIAGRAMS 

In order to determine age, metallicity, distance, mean Galac- 
tic reddening, and binary fraction we have compared the 
observational CMD with a library of artific ial populations 
Monte Carlo generated (see lTosi et al.ll2007r ). Different sets 
of stellar tracks[f] have been applied, fitting both primary 
age-sensitive descriptors, like the luminosity of the MS red- 
dest point ("red hook", RH), the red clump (RC) and the 
MS termination point (MSTP, evaluated at the maximum 
luminosity reached after the overall contraction, OvC, and 
before the runaway to the red), and secondary CMD fea- 
tures (sensitive to a broad range of parameters) like the RH 
colour, the luminosity at the base of the red giant branch 
(RGB), the RGB colour and inclination, the RC colour. The 
most valuable age indicator is in principle the TO point, 
evaluated at the bluest point after the overall contraction; 
however, for young and intermediate age clusters this phase 
is poorly populated. More in general, colour constraints are 
less reliable than luminosity constraints, since the former 
are much more influenced by theoretical uncertainties like 
colour transformations and the super-adiabatic convection. 

In order to make a meaningful comparison, all synthetic 
CMDs are combined with stars picked from an equal area of 

3 The same used in our past works, to maintain homogene - 
ity: the old Pa dova jBressan et ah! Il993l ; iFagotto etahl Il994h . 
the FRANEC llDominguez et alj I1999T). and t he FST ones 
dVentura et~al]|l998h . See lBragaglia fc Tosil j2006ll for a descrip- 
tion of their main properties and a detailed justification of their 
use even if newer tracks have appeared in the meantime. 
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Table 3. New clusters structural parameters; numbers in parenthesis are the errors on RA and Dec. 



Cluster 


RA Dec 


RA Dec 


r c 


rh 


n 


c 




present paper 


WEBDA 










Be 23 


06 33 15 (1.24") 20 31 57 (1.98") 


06 33 30 20 33 00 


55" 


178" 


1458" 


1.4 


Be 31 


06 57 37 (2.84") 08 18 20 (2.91") 


06 57 36 08 16 00 


50" 


171" 


1540" 


1.5 


King 8 


05 49 18 (3.21") 33 37 50 (2.22") 


05 49 24 33 38 00 


60" 


89" 


325" 


0.7 



the adjacent field which is located about 9' away from the 
cluster centre. For each cluster, the synthetic CMD is pop- 
ulated until the total number of stars (synthetic plus field) 
equals the observed number of cluster stars brighter than 
V — 20 and corrected for photometric errors and incom- 
pleteness, as derived from the artificial star tests. 

As a first step, we seek cluster parameters that are com- 
mon to all solutions. Differential reddening and fraction of 
binaries are the typical cases, because they are almost in- 
dependent of the specific set of stellar tracks. When these 
quantities are fixed (essentially matching the MS width), 
the magnitude difference between RH, MSTP and RC is 
effectively used to co nstrain the age (see e.g., Fig. 8 in 
ICastellani et al.ll2003l ). Once the best fitting age is estab- 
lished, all remaining parameters (mean Galactic reddening, 
metallicity, and distance) are varied till the colour of the 
MS and the magnitudes of RH, MSTP, RC are matched. A 
final test of the fit quality considers the RH, RGB and RC 
colours as well as the luminosity at the base of the RGB. 

Critical points of this process are the MS and the RC 
fitting. The former issue involves the upper and the lower 
ends of the MS. The RH morphology can vary from a verti- 
cal orientation to a very hooked morphology with the change 
of age, metallicity, overshooting, microscopic physics, etc., 
while the inclination of the lower MS follows the equation of 
state and the adopted colour transformations. In this con- 
text, a combination of parameters is defined "acceptable" 
when the corresponding synthetic CMD fits "most" of the 
visible MS shape. On the other hand, the RC fitting can 
be hampered by both theoretical and statistical uncertain- 
ties. The RC morphology and luminosity strongly depend 
on fundamental parameters like the metallicity, the age and 
the helium content, as well as by more subtle physical inputs 
like the efficiency of core overshooting an d the amount of 
mass loss in the pre He burning phase (see ICastellani et al.l 
2000 for an in-depth discussion). Moreover the number of 
RC stars is rather small in all three clusters, increasing the 
probability of confusion with RGB and field stars and not 
allowing very precise constraints on the derived age. 

In the following, we present the three clusters from the 
oldest to the youngest one. Note that we did not make use of 
the existing information on RVs for Be 31 because they are 
inconclusive regarding the me mbership (see Sect. 1 and the 
discussion in lFriel et al-l feoiO) and for Be 23 because they 
are limited to two stars. 



tain a significant number of cluster members, we restrict our 
sample to the stars within a radius of 2.5' from the cluster 
centre (that is 2.5 times r c , slightly more than ru, and well 
within r t , see Fig. [5] and Table [3}. The corresponding CMD 
is used for the determination of the cluster parameters. The 
most important features are emphasised in the left panel of 
Fig. [6] (while the right panel shows the control field CMD). 
First of all, the OvC gap (between V = 17.45 and V = 17.6) 
and the base of the RGB (around V = 17.65) are clearly vis- 
ible. The region of the RC is not well populated, but it is at 
least identifiable with the mild, but in our view significant, 
excesfl of stars in the magnitude range 15 < V < 15.6 and 
colour B - V » 1.1 - 1.2. 

Finally, an excess of stars running parallel and red- ward 
to the MS indicates that a fraction of stars is in binary 
systems. 

A remarkable number of stars is located blue-ward 
(B — V <0.4) and mainly above (15 < V < 17.5) the MSTP 
(see Fig. [6} , which cannot be explained by field stars (the 
control field does n ot contain any such blue obje ct). As al- 
ready suggested by lAhumada fc Lapassetl (|2007l ) these ob- 
jects are likely blue stragglers. 

As already found several times in the past, the MS is 
broader than we would have expected from the photomet- 
ric errors; what is surprising is that such an effect is larger 
around the RH. This might be a result of the combined ef- 
fect of differential reddening and RH morphology. On a more 
speculative side, keeping in mind the strong dependence on 
age of the RH morphology, it is also conceivable that the 
RH colour width might represent a signature of a prolonged 
star formation. 

To solve the problem of the MS scatter, one has to con- 
clude that either a differential reddening of about E(B — 
V) — 0.05 is present or a large fraction of binaries hosts 
very low mass companions (the photometric properties of 
a binary system are very close to the one of the primary 
alone, if the primary is below 2Mq and the companion's 
m ass is less than half the primary's mass; see e.g. Fig. 1 
in lHurlev fc Toutlll998l ). or both. Given the patchy nature 
of the obscuring material in the Galactic thin disc we are 
more inclined to the reddening hypothesis. In addition to 
the mean Galactic reddening, in the following we will adopt 
a differential reddening of about 0.05 mag unless otherwise 
stated. 

A rough estimate of the binary fraction can be obtained 



4.1 Berkeley 31 

In this group of OCs, Be 31 is the richest in term of mem- 
bers. The contamination of field stars is evident, especially 
above the sub-giant branch and blue-ward of the RGB. In 
order to reduce the number of field interlopers, but still re- 



4 Such an excess is found with respect to other equal-area CMD 
randomly chosen from the control field. 

5 Here the total reddening is considered a sum of a mean Galac- 
tic component and a random component ("differential") which 
includes the foreground fluctuation and a possible internal red- 
dening. 
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Figure 6. Left panel: CMD of stars within 2.5' of the centre of Be 31. We indicate the luminosity level of the RH, the MSTP, the RC 
and the base of the RGB (BRGB). The dotted and the dot-dashed boxes are, respectively, used to estimate the fraction of single and 
binary stars. Right panel: CMD of field stars (as evaluated from an equal area adjacent region). 



by evaluating the number of stars in two appropriate CMD 
boxes, one along the MS (see the dotted line in Fig. [6} and 
one red-ward of the MS (chosen to cover the binary sequence; 
see the dot-dashed line). We apply magnitude cuts and keep 
only stars below the RH and brighter than V — 22 to avoid 
evolved stars where the binary evolution can be quite com- 
plicated. To remove the field contamination we have statisti- 
cally subtracted the contribution of field stars using several 
equal area regions randomly chosen from the control field. 
The net result is a binary fraction between 22% and 26%. 
Strictly speaking these estimates may be slightly underesti- 
mated with respect to the real fraction, since low mass ratio 
binaries, whose properties are close to those of single stars, 
might be partially missed. However, a mean fraction of 24% 
appears to be a reasonable "ansatz" and will be assumed for 
all the simulations. 

Keeping the binary fraction and the differential redden- 
ing fixed, we have searched the best combination of param- 
eters for each set of tracks and metallicity. Table [4] sum- 
marises the results and the residual discrepancies. Synthetic 



RGB, RC, and lower MS colours together with the magni- 
tude of the RGB base are labelled as "TB" , "TR" , "TF" , 
or "OK" when they are too blue, too red, too faint, or in 
agreement with respect to the data. 

Concerning the age, the overall interval of confidence 
is between 2.3 and 2.9 Gyr. As expected, to reproduce the 
observed colours, models of progressively lower metallicity 
require a larger mean reddening and a larger distance mod- 
ulus, because they are intrinsically bluer and brighter. How- 
ever, the quality of the fit subtly depends on the adopted 
metallicity. 

The FST models with Z = 0.006 and Z = 0.01 fit rea- 
sonable well the RH, MSTP and RC luminosity levels, the 
colour of the RGB and the magnitude of the RGB base but 
they predict a "too hooked" RH (that doesn't mean the 
colour is wrong, but rather that the MS shape is too curved 
before the RH point), a synthetic MS (below the 19th mag- 
nitude) and a RC bluer than actually observed. In term of 
cluster parameters, the former metallicity would suggest a 
cluster age of 2.6 Gyr, E(B - V) = 0.165 ± 0.025 and a 
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Table 4. Be 31 best fit parameters (age, mean Galactic reddening, differential reddening, and distance modulus) for various sets of 
models and metallicities. Columns from 7 to 10 indicate how the models reproduce the following CMD features: RGB colour, magnitude 
at the base of the RGB, RC colour, and lower MS colour. Additionally, the 10th column reports the magnitude below which the synthetic 
MS diverges from the data MS. See text for definition of the flags TR, TB, TF, OK. 



Set 


Z 


Age (Gyr) 


E(B - V) M 


E(B-V) D 


(m - M) 


RGB (C) 


RGB base (M) 


RC (C) 


Low MS (C) 


FST (»7 = 0.2) 


0.02 


2.7 


0.020 


±0.00 


14.90 


TR 


OK 


TB 


TB(V>19) 


FST \n = 0.2) 


0.01 


2.8 


0.095 


±0.025 


14.60 


OK 


OK 


TB 


TB(V>19) 


FST (ij = 0.2) 


0.006 


2.6 


0.165 


±0.025 


14.60 


OK 


OK 


TB 


TB(V>19) 


Padova 


0.02 


2.8 


0.00 


±0.00 


14.75 


TR 


TF 


OK 


TB(V>20.5) 


Padova 


0.008 


2.6 


0.135 


±0.025 


14.45 


TR 


OK 


TB 


TB(V>20.5) 


Padova 


0.004 


2.9 


0.185 


±0.025 


14.40 


OK 


OK 


TB 


TB(V>20.5) 


FRANEC 


0.02 


2.5 


0.025 


±0.025 


14.80 


TB 


TF 


TR 


OK(V>21) 


FRANEC 


0.01 


2.3 


0.145 


±0.025 


14.60 


TR 


OK 


OK 


OK(V>21) 


FRANEC 


0.006 


2.5 


0.175 


±0.025 


14.46 


OK 


OK 


TB 


OK(V>21) 



distance modulus (m — M) a = 14.6, the latter a cluster age 
of 2.8 Gyr, E(B - V) = 0.095 ± 0.025 and a distance mod- 
ulus (m — M) = 14.6. Raising the metallicity to Z = 0.02 
worsens significantly the quality of the fit. In fact, the mag- 
nitude difference between RC and MSTP requires an age 
of 2.7 Gyr and any attempt to reconcile the predicted and 
observed RH colour and magnitude runs into two additional 
problems: the reddening should be l owered to an untena ble 
E(B - V) « 0.02 (compared to the ISchlegef et al. I ll998l es- 
timate of 0.145) and the synthetic RGB would be too red. 

With the Padova stellar tracks, we find that the lower 
metallicity, Z — 0.004, provides the best match, although 
suffering from the same pathologies noticed with the Z = 
0.006 FST set, i.e., too blue synthetic RC and lowest MS. 
In this case the estimated age raises to 2.9 Gyr, and we find 
a mean reddening E(B — V) = 0.185 ± 0.025 and a distance 
modulus (m - M) a = 14.4. The metallicity Z = 0.008 does 
not improve the fit and makes the match worse, producing 
a too red RGB. Accepting these discrepancies we get an 
estimated age of 2.6 Gyr, E(B - V) = 0.135 ± 0.025 and 
distance modulus (m — M) — 14.45. Finally, the Padova 
Z — 0.02 set is not acceptable for the same reason as the 
solar FST: matching the MSTP, RC and the RH magnitudes 
would require a null reddening, contrary to ISchlegel et al.l 
i| 19981 ) values. 

The FRANEC models are the only ones able to repro- 
duce well the faint end of the MS (down to V « 20 - 20.5). 
Using the metallicity Z — 0.006 we can reproduce all impor- 
tant features apart for a mismatch in the colour of the RC 
that is too blue. In this case the best combination of param- 
eters is: cluster age 2.5 Gyr, E(B - V) = 0.175 ± 0.025, 
and distance modulus (m — M) — 14.46. The metallic- 
ity Z = 0.01 reduces the RC colour discrepancy, but to 
the cost of a synthetic RGB that is too red by as much 
as ~ 0.1 mag and under-luminous at the base by as much 
as ~ 0.25. Using Z — 0.01 the best set of parameters is: 
age of 2.3 Gyr, E(B - V) = 0.145 ± 0.025, and distance 
modulus (m — M) — 14.6. A very similar result can be 
achieved by using a metallicity Z = 0.02, provided that all 
the main parameters be readjusted: age 2.5 Gyr, differential 
reddening E(B ~ V) = 0.025 ± 0.025 and distance modulus 
(m - M) = 14.8. 

Figure [7] summarises the best fitting CMD for each set 
of tracks, compared with the observational CMD (top left 
panel) . 



Trying to understand the results in terms of different 
prescriptions of the input physics, a first interesting con- 
clusion concerns the efficiency of the core overshooting. As 
expected, the RH length is severely affected by this macro- 
scopic effect as a consequence of the longer core hydrogen 
burning. Models with overshooting predict ages higher (2.6- 
2.9 Gyr) than those inferred with canonical models (2.3-2.5 
Gyr). Our analysis suggests that the observed RH morphol- 
ogy is better delineated by the non-overshoot tracks, as rep- 
resented by the FRANEC models, while the Padova and 
FST solutions, which use different amounts of overshooting, 
are too hooked. On the other hand, the low MS is unaf- 
fected by overshooting, being these stars characterised by 
radiative cores. In this case, the good match offered by low 
mass FRANEC models reflects, rather, a combined effect of 
the equation of state and of atmosphere model. 

On the basis of this discussion, the FRANEC set should 
be preferred, since it is in agreement with both the RH mor- 
phology and with the low MS location. This restricts the set 
of possible ages to the range from 2.3 to 2.5 Gyr. Concern- 
ing the metallicity, a solar abundance is ruled out by a clear 
inability to reproduce the RGB colour. This result is in rea- 
sonable agreement with what has been found in literature, 
both by photometric and spectroscopic methods (see Sect. 
!)• 

With these ranges of age and metallicity, the mean 
Galactic reddening is bet ween 0.145 and 0.175 mag, only 
slightly higher than the ISchlegel et al.l |l998) estimate 
E(B — V) = 0.145 (which is an upper limit, being the asymp- 
totic reddening in the cluster direction, but which is also 
more uncertain at low Galactic latitudes like those of our 
clusters), while the best estimate for the distance modulus 
is between 14.46 and 14.6. 

A final comment concerns the choice of adopting a frac- 
tion of binaries of 24%: surprisingly, a visual inspection of 
Fig.[7]seems to suggest that this number is slightly too high. 
However, this is clearly an artifact of the minimum mass 
(0.6 Mq) in our stellar models, whose effect is to bias the 
binary distribution toward the high mass ratio. Although 
this effect is negligible near the TO, it grows near the lower 
end of the MS. 
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Figure 7. Top left panel: CMD of stars within 2.5' of the centre of Be 31. The other panels, clockwise from this, show the best fitting 
synthetic CMDs for the following parameters: FRANEC Z = 0.01, age 2.3 Gyr, E(B - V) = 0.145 ± 0.025 and (m - M) = 14.6; Padova 
Z = 0.004, age 2.9 Gyr, E(B - V) = 0.185 ± 0.025 and (m - M) = 14.4; FST Z = 0.01, age 2.8 Gyr, E(B - V) = 0.095 ± 0.025 and 
(m — M)o = 14.6. The adopted percentage of binaries (with random mass ratio) is always 24%. 



4.2 Berkeley 23 



The strong contamination and the apparent small number of 
cluster stars make the CMD of Be 23 rather poorly defined. 
In the left panel of Fig. [S] we show the CMD restricted to 
the inner 3' from the cluster centre (that is about three 
r c and slightly less than r^, see Fig. [5] and Table [3j) . This 
partially cleans the diagram, revealing some signatures of 
the cluster: 1) a broad and irregularly clumped MS, clearly 
distinguishable at least down to V = 22; 2) a RH positioned 
near V = 16.2 (evidenced in blue); 3) a MSTP near V = 15.7 
(in red); 4) a mild excess of stars which is compatible with 
an apparently elongated RC, extending between V = 15.1 
and V = 15.5 (green). The RGB, not sufficiently populated 
to be distinguishable from the field, is instead difficult to 
recognise. This circumstance, together with the wide OvC 



gap and the extended RC are indications that the cluster is 
younger than 2 Gyr. 

Before applying the synthetic CMD approach, we have 
evaluated the amount of differential reddening and the frac- 
tion of binaries. The MS presents a substantial intrinsic scat- 
ter, larger than expected from photometric errors, suggest- 
ing a differential reddening of at least 0.05 mag. The percent- 
age of binaries is not easy to obtain, given the strong field 
contamination. Using the same strategy applied to Be 31, we 
have found a range of values between 20% and 30%. Such 
large uncertainty should be taken as a consequence of the 
spatial fluctuations across the control field. In the following 
we adopt a differential reddening E(B — V) = 0.05 and a 
mean fraction of binaries of 25%. 

In order to put limits on the cluster age and metal- 
licity, the 3' CMD is compared with our synthetic CMDs. 
Generally speaking, all examined models do not satisfacto- 
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Figure 8. Left panel: CMD of stars within 3' of the centre of Be 23. Labels indicate the luminosity level of the RC, the MSTP and the 
RH. Right panel: control field CMD (taken from an equal area adjacent region). 



rily reproduce two evident CMD features. First of all, the 
observed RH morphology is more vertical than predicted by 
the models. Furthermore, we cannot effectively reproduce 
the number counts at the MSTP and at the RC, which are 
always under-predicted by all models. 

Despite these discrepancies, some conclusions can be 
drawn using the RH, MSTP and RC luminosities as well as 
the MS colour. Table \S\ summarises the best combination 
of parameters and the residual discrepancies for each set of 
tracks. 

All the examined models can be reconciled with the lu- 
minosity constraints, whereas only the FST and the Padova 
tracks at solar metallicity can fit the average RH colour and 
the RC colour. In both cases the cluster age is found to 
be around 1.3 Gyr. At lower metallicities the synthetic RH 
colours become systematically too blue. This drawback is 
particularly severe for the FRANEC tracks. The reason for 
this is the younger ages required for the FRANEC set to fit 
the luminosity constraints (RH, MSTP, RC), since they do 
not consider overshooting. However, the FRANEC models 
work better than any other to reproduce the colour of the 



lower MS. In particular, these tracks are able to match the 
MS down to V « 21 while all other models diverge to the 
blue below V « 20. 

Figure shows the best fitting CMD for each set of 
tracks and the corresponding parameters. We have decided 
to give preference to the RH region, which provides infor- 
mation on the age, rather than the low MS fit, hence all 
FRANEC models have been rejected. With this assump- 
tion, the age for Berkeley 23 is well constrained between 1.1 
and 1.3 Gyr, while the favoured metallicity is solar. Con- 
cerning reddening and distance, the range of variation for 
the former is 0.225-0.425 and for the latter is between 13.6 
and 14.00. Furthermore, it is noteworthy that the reddening 
required for the solar models ( E(B - V) = . 225 ± 0.025) 
is too low compared with the ISchlegel et ail |l998) value 
(E(B - V) = 0.35), while that implied by the half solar 
models is perfectly consistent, and the reddening required by 
even metal poorer models (Z=0.006 and Z=0.004) is barely 
acceptable or too high. Combining these arguments we sug- 
gest that the metallicity of Be 23 is probably between solar 
(as suggested by the CMD fitting) and half solar. 
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Table 5. Be 23 best fit parameters (see previous table). 



Set 


Z 


Age (Gyr) 


E(B - V) u 


E{B - V) u 


(m - M) 


RH(C) 


RC (C) 


Low MS (C) 


FST (77 = 0.2) 


0.02 


1.3 


0.225 


±0.025 


14.00 


OK 


OK 


TB (V>20) 


FST (jj = 0.2) 


0.01 


1.2 


0.325 


±0.025 


13.80 


TB 


OK 


TB (V>20) 


FST (tj = 0.2) 


0.006 


1.1 


0.385 


±0.025 


13.75 


TB 


OK 


TB (V>20) 


Padova 


0.02 


1.3 


0.225 


±0.025 


14.00 


OK 


TR 


TB (V>20) 


Padova 


0.008 


1.3 


0.325 


±0.025 


13.75 


TB 


OK 


TB (V>20) 


Padova 


0.004 


1.2 


0.425 


±0.025 


13.6 


TB 


OK 


TB (V>20) 


FRANEC 


0.02 


1.0 


0.245 


±0.025 


13.95 


TB 


OK 


TB (V>21) 


FRANEC 


0.01 


0.9 


0.375 


±0.025 


13.75 


TB 


OK 


TB (V>21) 


FRANEC 


0.006 


0.8 


0.445 


±0.025 


13.75 


TB 
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Figure 9. Top left panel: CMD of stars within 3' of the centre of Be 23. The other panels, clockwise from this, show the best fitting 
synthetic CMDs for the following parameters: FRANEC Z = 0.01, age 0.9 Gyr, E(B-V) = 0.375 ±0.025 and (m-M) = 13.75; Padova 
Z = 0.008, age 1.3 Gyr, E(B - V) = 0.325 ± 0.025 and (m - M) = 13.75; FST Z = 0.01, age 1.2 Gyr, E(B - V) = 0.325 ± 0.025 and 
(m — A'/)q = 13.8. The adopted percentage of binaries (with random mass ratio) is always 25%. 
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As already noted in Sect. 1, there is discrepancy be- 
tween different analyses of this cluster, and our st udy makes 
no exception. We disagree with lAnn et alJ ({2002) , who find 
the cluster to be much younger and metal-richer, yet more 
reddened, even if these authors used MSTP and RH lu- 
minosities very close to the values adopted here (see their 
Fig. 3). However, their best fit model seems to m i ss the RC 
position. We disagree also with lHasegawa et al.l (|2004T ) on 
the age and metallicity (their best solution is much older 
and metal-poorer) even if reddening and distance modulus 
are in very good agreement. The age and metallicity dif- 
ferences arise mainly from a different interpretation of the 
MSTP position: by looking at their Figure 3 it is evident 
that these authors have identified as the MSTP a point that 
is half a magnitude fainter than ours. 

4.3 King 8 

King 8 exhibits a very scattered CMD suggesting the 
strongest differential reddening among the three clusters. 
This is clearly evident from the left panel of Fig. [10] where 
the CMD of stars within 3' from the cluster centre (that is, 
three times r c , 2.5 rn and within r t , see Fig. [5] and Table 
3) is compared with the control field CMD (right panel). 
In these conditions few CMD features are available for the 
fit: 1) a MSTP around V = 16.0 (indicated with a red line 
in Fig. [10}; 2) a wide OvC gap; 3) a RC around V = 15.5 
(indicated with a green box). As for Be 23, there is no ev- 
idence of a well formed RGB. More importantly, the large 
MS spread prevents a safe identification of the RH, which 
is tentatively identified as the mildly hooked feature around 
V — 16.3 (indicated with a blue line). 

Given the wide (0.2 mag) and rather uniformly filled 
MS, the current data do not allow us to reach a firm con- 
clusion about the binary fraction. Hence, in the following 
simulations we have assumed a fixed fraction of 30% (with 
random mass ratio). On the other hand, a differential red- 
dening E(B — V) of at least 0.2 mag represents the only vi- 
able solution to reproduce the observed MS spread. Keeping 
fixed these parameters, we have investigated the possibility 
to fit simultaneously the MSTP, the RH and the RC lumi- 
nosities by adjusting the age, the mean Galactic reddening 
and the distance modulus. The best fitting results for the 
various sets and metallicities are displayed in Table [6] with 
the last two columns describing how close the models match 
the RC and the lower MS colours. 

In contrast with what is found for Be 23 and Be 31, 
all explored models match very well the lower MS (at least 
up to V — 21). On the other hand, as a consequence of 
the high differential reddening, all synthetic CMDs show a 
rather elongated RC while the observed RC is round and 
concentrated. Fig. [TT] summarises the best fitting CMD for 
each set of tracks compared with the data CMD (top left 
panel) . 

With the FST models, both the solar and the half- 
solar metallicities fail to reproduce the RC colour (bluer 
than observed) and the number of stars around the MSTP 
phase (systematically underestimated). A much better re- 
sult is achieved using the Z — 0.006 models, whose synthetic 
CMDs account for the RC colour, although still predicting 
too few MSTP stars. In terms of age, the Z — 0.006 case 
corresponds to about 1.1 Gyr, while the best fit distance is 



about (m — M)o = 12.95. Furthermore, the high reddening 
implied by the metallicity Z — 0.006 is in good agreement 
with the literature v alue (E(B — V) — 0.8) based on the 
ISchlegel et ail (|l998ft maps, while the reddening solutions 
E(B-V) = 0.55±0.10 and E(B-V) = 0.65±0.10, required 
by the metal richer tracks, are too low or barely acceptable, 
respectively. 

In contrast with the FST result, among the Padova 
models only those at solar metallicity predict the correct 
RC colour, whereas the synthetic RC is too blue (by about 
0.1 mag) at Z = 0.008 and Z = 0.004. On the other hand, no 
significant difference is found in the best fit age with vary- 
ing metallicity: with Z = 0.02 the best CMD is for 1.3 Gyr, 
while with Z — 0.004 it is 1.1 Gyr, results which are almost 
identical to the FST findings. Focusing on the MSTP, both 
Padova and FST tracks share the same mismatch (which is 
to be ascribed to the overshooting): the predicted number 
of MSTP stars is always underestima ted. As found f or the 
FST models, the comparison with the ISchlegel et ail (|l998T l 
maps reddening favours the lowest metallicity (Z — 0.006), 
in contrast with our photometric estimate. 

As far as the FRANEC models are concerned, we don't 
find a clear indication to prefer one metallicity over the 
others. Both the observed MSTP counts and the lower MS 
colours are fairly well reproduced (the latter down to V w 
22), while the mismatch here concerns the predicted number 
of RC stars (whose colours are in good agreement with the 
observations), which outnumber the observed counts by a 
factor of 3. The resulting ages (0.8 - 0.9 Gyr) are younger 
than those obtained with the Padova and FST models, as a 
direct consequence of the lack of overshooting. Once again, 
only the sub-solar models (Z = 0.01 and Z = 0.006) imply 
a reddening value compatible with the lSchlegel et all (|l998T ) 
estimate. 

Also for this cluster the comparison to previous results 
is not simple. We seem to agree, apart from the Padova so- 
lution, on a rather low met a l cont ent. The apparent agree- 
ment with IChristianl (|l98ll , Il984l ) is not significant, given 
the differences betwe en the two photometries (see Fig. 3}. 
IKoposov et al.l |2008T ) adopted solar isochrones for all their 
clusters and this could at least partly explain the differences, 
although the dependence on metallicity should be smaller in 
the near IR bands. 



5 CONCLUSIONS 

5.1 Galactic considerations 

Be 31 and Be 23 are rather far away from both the Galactic 
centre and the Galactic plane. Berkeley 31 has a Galactocen- 
tric distanc^fl (Rgc) between 15.4 kpc and 15.9 kpc and a 
height above the Galactic plane between 693 pc and 739 pc, 
while Berkeley 23 has a Rgc between 13.1 kpc and 14.2 kpc 
and a height between 494 pc and 594 pc. Such scale heights 
could be the result of a formation in situ as a part of the 
Thick Disc population, as also suggested by the old age of 
Be 31 (well over 2 Gyr), or due to a disc flaring which grad- 
ually increases the Galactic scale-height of the Thin Disc. 
On the other hand, the shorter scale height (between 196 pc 

6 Assuming a Galactocentric distance of the Sun of 8 kpc. 
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Figure 10. Left panel: CMD of stars within 3' of the centre of King 8. Right panel: control field CMD (as evaluated from an equal area 
adjacent region). 



Table 6. King 8 best fit parameters. Columns 7 and 8 indicate the ability of the models to reproduce the RH colour and the lower MS 
colour. Additionally, the 8th column reports the magnitude below which the synthetic MS diverges. 
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and 236 pc), its Rgc between 11.6 kpc and 12.3 kpc and 
the younger (probably less than 1 Gyr) age make King 8 a 
likely member of the Thin Disc. 

Such findings help to interpret the difference in the field 
population between the direction of Be 23 and Be 31 and the 
direction of King 8 visible in the lower panels of Fig. [1] In the 



former, a structure to the blue of the clusters' main sequence 
is clearly visible starting from V = 19. Such a contamination 
is a signature of a much farther and probably older popula- 
tion along the line of sight. The most l ikely candidate is th e 
above mentioned Thick Disc (see e.g. ICignoni et alj [20081. 
whose scale height is at least three times longer than the 
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Figure 11. Top left panel: CMD of stars within 3' of the centre of King 8. The other panels, clockwise from this, show the best fitting 
synthetic CMDs for the following parameters: FRANEC Z = 0.01, age 0.9 Gyr, E(B - V) = 0.67 ± 0.10 and (m - A/) = 13.0; Padova 
Z = 0.02, age 1.3 Gyr, E(B - V) = 0.52 ± 0.10 and (m - M) = 13.2; FST Z = 0.006, age 1.1 Gyr, E(B - V) = 0.75 ± 0.10 and 
(m — M)o = 12.95. The adopted percentage of binaries (with random mass ratio) is always 30%. 



corresponding Thin Disc. Although less probable, it may 
also represents a tida l debris of a disrup ted galaxy (see e.g. 
the Monoceros ring, iMartin et alJ l2004h or a signature of 
the Galactic Warp and Flare ( Momanv et alj|200a ). On the 
other hand, the lack of this feature in the direction of King 8 
is mostly due to the closer proximity to the Galactic plane 
(3 degrees against ~ 5 degrees of Be 23 and Be 31). At 
this latitude the line of sight goes deeper into the Thin Disc 
(for about 6 kpc) before crossing a pure Thick Disc sample, 
making the Thin Disc stars the dominant contaminant. Fur- 
thermore, the Galactic reddening increases with distance; 
this explains the elongated shape of the M-dwarfs sequence 
(around B - V ~ 1.9) exhibited by the King 8 CMD field, a 
circumstance that is not observed either in Be 23 or in Be 31 
(where the M-dwarf sequence is rather vertical). 



5.2 Summary 

Our study of open clusters is aimed at better understanding 
the chemical and structural properties of the Galaxy. In this 
paper we have studied three distant open clusters toward 
the Galactic anti-centre using deep LBT photometry. The 
CMDs resulting from these data are more precise and more 
than three mag deeper than the ones found in literature, ex- 
tending to more than six mag below the MS TOs. The syn- 
thetic CMD technique has been used to derive a confidence 
interval for age, metallicity, binary fraction, reddening, and 
distance for each cluster. To remove the model dependence 
of the results three different sets of stellar tracks (Padova, 
FST, FRANEC) have been adopted. From the comparison 
between models and data we have drawn the following con- 
clusions: 



• Be 31 is located at about 15-16 kpc from the Galactic 
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centre and about 700 pc above the Galactic plane. The re- 
sulting age varies from 2.3 Gyr to 2.9 Gyr, depending on the 
adopted stellar model, with better fits for ages between 2.3 
Gyr and 2.5 Gyr. Concerning the metallicity we obtained 
a good match using the sets of models with metal content 
lower than solar. The mean Galactic reddening E(B — V) is 
between 0.095 and 0.185 mag and the fraction of binaries is 
between 22% and 26%. 

• Be 23 is at about 13-14 kpc from the Galactic centre 
and 500-600 pc above the plane. The age is well constrained 
between 1.1 and 1.3 Gyr, while the best fitting metallic- 
ity is more uncertain: the CMD fit suggests a solar value, 
while only a half-so lar metallicity is compatible with the 
ISchlegel et al.l (| 19981 ) estimate for the Be 23's direction. The 
mean Galactic reddening E(B — V) is between 0.225 and 
0.425 and the binary fraction is between 20% and 30%. 

• King 8 is at about 12 kpc from the Galactic centre and 
about 200 pc above the plane. The strong differential red- 
dening up to ± 0.1 mag broadens the colour extension of the 
MS, hindering a precise estimate of the cluster parameters. 
In this condition, the best fitting age is 0.8-1.3 Gyr, while the 
entire range of metallicities (0.004-0.02) is consistent with 
the data. On the other hand, only sub-solar models lead 
to reddening estimates (ranging fro m E(B — V) = 0.55 t o 
E(B ~V) = 0.88) in agreement with ISchlegel et all l|l998Tl . 
A 30% of binaries is compatible with the observed CMD. 

With these results, Be 31 and Be 23 are candidate Thick 
Disc open clusters, while King 8 is a more classical Thin 
Disc member. Only further direct spectroscopic and kine- 
matic searches will allow to shed light on the Galactic origin 
of these clusters. They will also be very important for any 
study of the radial distribution of metallicity (and chemi- 
cal abundances) in the Galactic disc, a crucial ingredient in 
chemical evolution models. 

Our choice to exploit the RC position even if it is barely 
visible in the three clusters deserves a final comment. For 
Be 31, which shows the least defined RC, many other rel- 
evant and well defined CMD features (RH, MSTP, RGB, 
sharp MS) are available. Hence, even without using the RC 
we would get t he sa me parameter estimates. Moreover, the 
ISchlegel et alj (1 19981 ) estimate of the reddening in this direc- 
tion can be safely used as an additional constraint. On the 
other hand, the CMDs of Be 23 and King 8 are both suffi- 
ciently contaminated and affected by reddening that the RC 
position is a necessary constraint to derive their ages. Re- 
laxing the assumption on the RC position greatly expands 
the range of possible ages and the only derivable result is 
that both clusters are certainly younger than 2 Gyr. Still, 
we believe that, although admittedly few, the stars in the 
CMD RC region are most likely cluster members and can 
therefore be used in the analysis. A definitive confirmation 
clearly needs proper motion or radial velocity studies. 

With the present study we may have reached the limit 
of what is possible to derive simply with photometric data. 
Next step will be provided by Gaia, the satellite due to fly in 
about two years, which will produce photometry for about 
10 9 Galactic objects, and most importantly, parallaxes and 
proper motions of unprecedented precision. The expected 
performances of Gaia will permit to obtain individual dis- 
tances of RC stars better than 10% up to 8-10 kpc from the 
Sun, i.e., including almost the entire family of known OCs. 



The precision on proper motions will be even better (the 
estimated error on proper motions is about half the one on 
parallax), so that membership for a large fraction of stars 
will be available also for distant clusters. This will be impor- 
tant especially in cases like Be 31, where the separation in 
radial velocity between field and cluster stars is very small 
(see Introduction). 
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